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THE  EFFECT  OF  THE  CHOICE  OF  RESPONSE  MATRIX 
ON  UNFOLDED  BONNER  SPHERE  SPECTRA 


INTRODUCTION 

In  1960,  Bramblett  et  al  [1]  suggested  that  neutron  spectra  could  be 
determined  by  using  a  thermal  neutron  detector  in  the  center  of  several 
moderating  spheres  of  various  sizes.  Ibey  used  a  small,  cylindrical,  4  mm 
high  X  4  mm  dia.,  europium- activated  lithium  iodide  scintillation  crystal  as 
the  thermal  neutron  detector  in  the  centers  of  polyethylene  spheres  having 
dicuneters  of  2,  3,  5,  8,  and  12  inches.  Thermal  neutrons  arriving  at  the 
center  of  the  moderator  interact  in  the  scintillator,  primarily  by  the 
®Li(n,a)^H  reaction,  producing  a  4.79  MeV  «  particle  which  is  stopped  in  the 
crystal.  The  scintillator  is  coupled  to  a  photomultiplier  and  the  output  is 
processed  to  give  a  measure  of  the  thermal  neutron  fluence  at  the  center  of 
the  moderator.  Such  spheres  are  usually  called  "Bonner  spheres”  after  the 
senior  author  of  the  original  paper. 


Bonner  spheres  can  serve  as  a  spectrometer  because  the  combined  effect  of 
moderation  and  absorption  of  the  incident  neutron  spectrum  causes  the  number 
of  thermal  neutrons  arriving  at  the  centers  of  the  spheres  to  vary  with 
energy  in  different  ways  for  the  different  size  moderators.  The  response  of 
a  set  of  spheres  may  be  written 

E 

f.  max 

A.  =  I  a.(E)*(E)dE  j  =  1,2,***M  (1) 
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A.  is  the  counting  rate  of  the  jth  detector. 


is  the  response  function  of  the  jth  detector  as  a  function 
of  energy 
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♦(E)  is  the  neutron  fluence  as  a  function  of  energy,  and 
M  is  the  total  nvimber  of  detectors. 

Equation  (1)  is  known  formally  as  a  Fredholm  integral  equation  of  the 
first  kind.  Such  equations  can  only  be  solved  if  aj(E)  is  an  analytical 
function,  which  is  not  the  case  for  Bonner  spheres.  In  practice,  Eq.  (1)  is 
replaced  by  a  set  of  linear  equations  by  dividing  the  energy  region  into 
several  regions  having  constant  detector  responses  and  fluences.  Equation  (1) 
becomes 


N 

Aj  =  ^  ^  a^j^${k)  j  =  1,2, •••M 

k=l 


(2) 


where 

is  the  response  of  the  jth  detector  to  neutrons  in  the  kth 
energy  interval,  and 

N  is  the  total  number  of  energy  intervals. 

The  set  of  detector  responses,  called  the  response  matrix,  has  not 

been  measured  with  sufficient  accuracy  because  of  difficulties  in  obtaining 
and  characterizing  suitable  neutron  calibration  fields.  Hence,  calculated 
response  matrices  are  usually  employed.  Several  such  matrices  have  been 
published,  with  three  appearing  in  the  past  few  years. 

The  purpose  of  our  study  was  to  determine  how  the  choice  of  response 
matrix  affects  the  neutron  spectrum  obtained  from  Bonner  sphere  data.  In 
particular,  we  were  interested  in  how  the  choice  of  matrix  affects  integral 
parameters  which  are  calculated  from  the  unfolded  spectrum.  These  parauneters 
include  dose,  dose  equivalent,  total  fluence,  quality  factor,  average  energy, 
and  the  response  of  personnel  monitoring  devices  and  radiation  survey 
instruments.  They  are  calculated  using  the  equation 


P  = 


N 

^  Ck$(k) 
k=l 


(3) 


«-• 


2 


where 


P  is  the  parameter  of  interest, 

$(k)  is  the  neutron  fluence  in  the  kth  energy  interval,  and 

C]j  is  the  conversion  factor  per  neutron  in  the  kth  energy  interval. 

THE  RESPONSE  MATRICES  USED  FOR  THIS  STUDY 

when  Bramblett  et  al  [1]  published  the  original  work  on  the  Bonner  sphere 
spectrometer,  they  reported  the  experimentally  determined  response  of  the 
detectors  for  neutrons  having  energies  between  0.05  and  15.1  MeV.  Measure¬ 
ments  were  also  made  for  thermal  neutrons,  and  the  response  of  the  various 
size  moderators  was  estimated  for  neutrons  having  energies  between  0.025  eV 
and  50  keV,  since  no  neutron  sources  were  then  available  having  energies  in 
this  region.  Using  this  data,  Burrus  [2]  reported  a  response  matrix  having 
9?  energy  bins  from  thermal  to  200  MeV.  This  matrix  is  usually  referred  to 
as  M60  [3] . 

Using  an  adjoint  transport  technique,  Hansen  and  Sandmeier  [4]  calculated 
the  response  of  the  original  set  of  Bonner  spheres.  These  calculated  responses, 
plus  those  for  moderators  having  diameters  of  10,  16,  18,  and  20  inches,  were 
reported  by  McGuire  [5].  A  52  energy-bin  version  of  these  calculations  was 
reported  by  O'Brien  et  al  [6].  This  52  bin  version  is  usually  referred  to  as 
M65  [3]. 

In  1973,  Sanna  [7]  calculated  a  31  group  set  of  responses  extending  from 
thermal  energies  to  400  MeV  using  the  adjoint  method  of  Hansen  and  Sandmeier. 

He  adjusted  the  ®Li(n,a)  cross  section  to  ensure  that  the  collision  probabil¬ 
ities  for  the  cylindrical  Li I  detector  would  be  retained  in  each  group  when 
it  was  represented  as  a  spherical  detector.  In  addition  to  making  calculations 
for  a  4  mm  X  4  mm  detector,  Sanna  made  calculations  for  8  mm  x  8  mm  and  12.7  ram  x 
12.7  mm  detectors  in  spheres  having  dieuneters  of  2,  3,  5,  8,  10,  12,  and  18 
inches.  He  also  investigated  the  effect  of  different  polyethylene  densities 
and  made  calculations  for  water  moderators  and  for  gold  foil  detectors.  For 
the  purposes  of  this  comparison,  we  call  the  response  matrix  for  the  4  mm  x 
4  mm  Lil  detector  in  polyethylene  spheres,  SAN4. 


3 


Recently,  three  new  response  matrices  have  been  caculcated.  Using  a 
Monte  Carlo  method,  response  matrices  were  calculated  at  the  Bhabha  Atomic 
Research  Center,  Bombay,  India  [8].  For  these  studies,  a  0.245  cm  radius  ^Lil 
scintillator,  100%  ®Li  enrichment,  and  a  density  of  4.061  gm/cm^,  was  taken  as 
being  equivalent  to  a  4  mm  x  4  mm  cylindrical  scintillator.  In  addition,  the 
effect  of  replacing  «iihe  Lil  scintillator  by  a  BF^  detector  was  studied.  We 
will  refer  to  the  matrix  for  the  4  mm  x  4  mm  detector  as  BARC. 


Based  on  experimental  data  and  calculations  available  in  the  literature, 
Zamborowski  [9]  showed  that  a  good  fit  to  this  data  could  be  obtained  by  as¬ 
suming  that  the  response  of  Bonner  spheres  of  various  sizes  was  given  by  the 
energy  response  function 


a  (E) 
J 


F(E) 


(E) 


(Inrj  -  M(E)) 


■] 


(4) 


S(E)rj  v/2? 


where 

Oj (E)  is  the  response  of  the  jth  sphere,  having  a  radius  in  cm  of  r, 
as  a  function  of  energy,  and 

F,  M,  and  S  are  parameters,  independent  of  sphere  size,  which  are 

adjusted  to  give  the  best  fit  to  the  calculated  and  experimental 
data. 

Values  of  F,  M,  and  S  were  tabulated  for  157  energies  between  0.4  eV  and 
15  MeV  for  a  4  mm  x  4  mm  Eil  detector.  The  response  matrix  obtained  from  these 
tabulated  values  and  Eq.  (4)  is  referred  to  as  LOGNM. 

The  most  recently  available  matrix  is  one  calculated  at  the  University 
of  Texas  at  Austin  by  llertel  and  Davidson  {10,11].  Using  newly  evaluated 
cross  sections  for  lithium,  carbon,  and  hydrogen,  the  method  of  Hansen  and 
Sandmeir  was  used  to  calculate  the  response  of  4  mm  x  4  mm  and  12.7  mm  x 
12.7  mm  Lil  detectors  and  for  moderators  having  diameters  of  2,  3,  5,  8,  10 
12,  and  18  inches.  Calculations  were  made  for  171  energies  between  thermal 
and  17.3  MeV  .  Calculations  were  also  made  for  the  bare  detectors  and  the 


detectors  in  2-,  3-,  and  5-inch  diameter  spheres  covered  with  0.76  cm 
cadmivim.  We  call  the  matrix  for  the  4  mm  x  4  mm  detector  UTA4. 

As  a  convenience  for  unfolding,  all  the  response  matrices  were  re-binned 
to  the  SANNA  energy  intervals.  For  M60  and  M65  this  was  done  by  making  a 
spline  fit  to  the  data  in  reference  [3].  The  neutrons  were  assumed  to  be  at 
the  log  center  of  the  SANNA  energy  intervals.  For  the  LOGNM  data,  the  re¬ 
sponses  were  calculated  for  each  of  the  157  energies  and  this  data  was  reduced 
to  the  SANNA  grid  using  log  interpolation.  Log  interpolation  was  also  done 
for  the  UTA4  data.  In  reducing  the  LOGNM  and  UTA4  data,  we  assumed  that  the 
original  energies  were  at  the  centers  of  log  energy  intervals.  Intervals  that 
fell  completely  within  a  SANNA  interval  were  converted,  weighted  by  the  ratio 
of  their  log  width,  to  the  SANNA  interval  log  width.  Intervals  which  crossed 
a  SANNA  interval  end-point  were  divided  using  a  log  scale  before  conversion  to 
the  two  SANNA  intervals.  No  numerical  data  was  available  for  the  BARC  matrix. 
Two  requests  to  the  authors  for  such  data  were  unanswered,  thus  the  data  for 
the  BARC  matrix  was  taken  from  the  graphs  in  reference  [8].  The  re-binned 
data  for  all  six  matrices  are  given  in  Tables  la-lf.  For  all  the  matrices, 
the  data  from  UTA4  is  used  to  calculate  the  effect  of  covering  the  spheres 
with  cadmium. 

THE  NEUTRON  SPECTRA 

The  neutron  spectra  and  the  associated  Bonner  sphere  data  used  to  test 
the  response  matrices  were  all  generated  at  the  Naval  Research  Laboratory  (NRL) 
during  the  course  of  establishing  the  response  of  the  Navy  albedo  neutron 
personnel  monitoring  badge  to  different  neutron  spectra.  Twenty  sets  of  data 
were  chosen,  covering  the  widest  range  of  neutron  energies  that  we  could  con¬ 
veniently  generate  using  Cf-252  and  AmBe  neutron  sources  moderated  by  Lucite, 
polyethylene,  steel,  and  cadmium  in  various  configurations.  This  data  was 
taken  using  a  set  of  Bonner  spheres  obtained  from  Ludlum  Measurements,  Inc., 
Sweetwater,  Texas.  The  detector  was  a  cylindrical  4  mm  x  4  mm  Lil  crystal. 

The  measured  density  of  the  polyethelene  was  0.96  g/cm^.  Total  counts  for 
each  detector  were  greater  than  10**  giving  statistical  errors  on  the  order  of 
±1%.  A  description  of  the  various  sources  is  given  in  Table  2,  while  the 
counts  for  the  2-,  3-,  5-,  8-,  10- ,  and  12-inch  detectors,  expressed  as  the 
percent  of  the  total  counts,  are  listed  in  Table  3.  The  spectra  have  been 


numbered  in  order  of  decreasing  average  neutron  energy.  Spectra  1,  4,  10,  and 
12  have  errors  considerably  greater  than  ±1%  because  they  were  separated  into 
source  and  room-return  components.  This  separation  was  done  using  the  method 
of  Eisenhauer  et  al  [12] .  Even  though  the  data  for  these  spectra  have  larger 
errors  because  of  this  process,  the  separation  allows  us  to  compare  the  experi¬ 
mental  data  to  relatively  well  known  calculated  data,  especially  for  Cf-252. 

Because  of  the  large  size  of  the  moderating  configurations,  it  is  not 
possible  to  separate  the  direct  and  room-return  components  of  these  spectra. 
The  effective  radius  of  the  experimental  room  is  included  as  an  aid  in  esti¬ 
mating  the  room-return  component.  The  room  having  an  effective  radius  of  3.4 
meters  measures  7.5  x  4.5  x  3.4  meters.  The  room  having  an  effective  radius 
of  5.0  meters  measures  12.2  x  6.1  x  4.6  meters. 

METHOD  OF  SPECTRUM  UNFOLDING 

The  neutron  spectra  were  unfolded  from  the  Bonner  sphere  data  using  the 
YOGI  neutron  unfolding  code  developed  at  NRL  by  Johnson  and  Gorbics  [13]. 

Since  data  for  the  bare  detector  were  not  available  for  all  matrices,  we  used 
only  the  data  for  the  2-,  3-,  5-,  8-,  10- ,  and  12-inch  detectors.  Only  the 
first  25  energy  intervals  were  used  for  the  spectrum  unfolding  since  data  were 
only  available  to  approximately  15  MeV  for  some  of  the  matrices.  Also,  from 
the  sources  used,  Cf-252  and  AmBe,  significant  numbers  of  neutrons  above  this 
energy  were  not  expected.  From  physical  considerations,  it  was  anticipated 
that  the  general  shape  of  the  spectra  could  be  described  by  a  three-component 
spectrum  made  up  of  a  high-energy  Maxwellian  peak,  an  intermediate  energy 
component  described  by  1/E^  (E  =  neutron  energy) ,  and  a  thermal  component. 
Therefore,  we  first  fit  the  data  using  the  MAXIET  algorithm  [14]  followed  by 
the  iterative  perturbation  method  of  YOGI .  We  did  not  smooth  the  spectra  to 
the  initial  MAXIET  spectrum.  However,  to  prevent  wild  oscillations  in  the 
spectra  due  to  errors  in  the  data  or  response  matrices,  we  stopped  the  unfold¬ 
ing  after  25  iterations,  using  an  8%  perturbation.  Further  unfolding  produced 
no  significant  improvement  in  the  fit  to  the  data. 

THE  CONVERSION  FACTORS 

The  conversion  factors  used  to  determine  integral  parameters  from  the 
unfolded  spectra  are  listed  in  Table  4.  The  response  of  the  Navy  albedo  badge 
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was  determined  by  making  the  best  fit  to  the  experimental  data  from  80  spectra 
generated  for  that  purpose.  The  response  of  a  completely  cadmium  covered 
albedo  badge,  often  called  a  "Hankins"  dosimeter,  was  approximated  from  the 
calculation  of  Allsmiller  and  Barish  [15].  The  response  of  the  Cr-39  detector 
was  taken  from  Benton  el  at  [16]  and  the  response  of  the  NTA  film  from  Oshino 
[17].  The  response  of  the  AN/PDR-70  was  supplied  by  Gordon  Riel  of  the  Naval 
Surface  Weapons  Center.  The  conversion  from  neutron  fluence  to  dose  was  taken 
from  the  tabulation  of  Sanna  [18],  while  the  conversion  to  dose  equivalent  was 
supplied  by  Charles  Eisenhauer  of  the  National  Bureau  of  Standards .  This  con¬ 
version  was  made  using  log-log  interpolation  [19]  of  the  data  in  ICRP  21  [20]. 

These  conversion  factors  may  not  necessarily  be  the  most  appropriate 
ones  available,  but  should  be  adequate  for  our  comparison  of  the  effect  of 
choice  of  response  matrix. 

RESULTS  AND  DISCUSSION 

General  Fit  to  the  Detector  Data 

One  measure  of  the  appropriateness  of  a  response  matrix  is  how  well 
the  Bonner  sphere  data  can  be  fit.  If  there  are  inconsistencies  in  the  re¬ 
sponse  matrix,  it  is  possible  to  have  data  which  cannot  be  fit  to  a  reasonable 
degree  of  accuracy  with  any  spectrum.  If  our  assumption  that  these  spectra 
can  be  fit  using  a  Maxwellian,  a  1/E^,  and  a  thermal  component  is  correct, 
then  the  fit  obtained  with  the  various  matrices  using  only  these  componenets 
is  an  indication  of  their  consistency.  The  root  mean  square  of  the  fit  to 
the  data  (%  average  error)  is  listed  in  Table  5  for  each  spectrum  unfolded 
using  each  matrix  and  the  MAXIET  algorithm.  Note  that  a  good  fit  is  obtained, 
especially  with  UTA4,  SAN4,  and  LOGNM.  There  appear  to  be  inconsistencies  in 
M60,  M65,  and,  especially,  BARC.  After  25  iterations  using  YOGI,  there  is 
considerable  improvement  to  the  fit  data,  the  average  average  error  being 
reduced  to  less  than  0.7%  for  SAN4  and  UTA4.  This  data  is  shown  in  Teible  6. 

In  order  to  illustrate  the  general  characteristics  of  the  fits  obtained  with 
each  matrix,  we  have  listed  the  percentage  difference  between  the  calculated 
response  of  each  detector  and  measured  response  in  Tables  7a-7f.  A  plus  sign 
indicates  that  the  measured  response  is  lower  than  the  calculated  response; 
but,  when  considering  the  matrices,  a  positive  difference  indicates  that  the 


calculated  response  of  the  matrix  is  too  high.  From  these  tables  we  can  draw 
the  following  conclusions : 


a.  For  M60,  the  responses  of  the  8"  and  12"  detectors  are  too  low, 
or  the  10"  response  is  too  high. 

b.  For  M65,  the  responses  of  the  3"  and  10"  detectors  are  too  high. 

c.  For  SAN4,  the  responses  of  the  2",  5",  and  12"  detectors  are 
slightly  low. 

d.  For  LOGNM,  the  responses  of  tiio  3",  8",  and  10"  detectors  are 
slightly  low,  or  the  other  three  are  too  high. 

e.  For  BARC,  the  responses  of  the  2",  5",  and  30"  are  too  low. 
Note,  particularly,  the  large  discrepancy  between  the  8",  10", 
and  12"  detectors. 

f.  For  U'rA4,  the  responses  of  the  2",  5",  and  12"  detectors  are 
slightly  low. 

These  conclusions  suggest  that  better  fits  to  the  data,  giving  more 
reasonably  shaped  spectra,  could  be  obtained  if  the  responses  of  selected 
detectors  were  raised  or  lowered  based  on  the  data  in  Tables  7a-7f.  We  did 
not  attempt  to  do  this  for  this  particular  study,  however. 

Agreement  witl^  Calculated  Spectra 

Spectra  4  and  12  were  chosen  for  comparison  with  the  calculations 
of  Grundl  and  Eisenhauor  [21]  and  Eisonhauer  et  al  [12].  In  Table  8  we  list 
the  calculated  counts  for  each  matrix  obtained  using  the  calculated  Cf-252 
spectrum  with  the  magnitude  of  the  spectrum  adjusted  to  minimize  the  percent 
difference  between  the  calculated  counts  and  the  experimental  counts.  Note 
the  remarkable  agreement,  especially  for  SAN4  and  UTA4 .  The  differences  indi¬ 
cate  that  M60,  M65,  hOGNM,  and  UTA4  will  give  spectra  having  lower  energies 
than  the  calculated  spectra.  Tiiere  is  no  discernable  trend  in  tiio  BARC  data. 
In  order  to  compare  the  experimental  spectra  with  tlie  calculated  spectrum, 
spectra  were  unfolded  using  each  matrix  witli  the  calculated  spectrum  as  a 
starting  spectrum  for  YOf;i .  'I’iie  unfolded  spectra  wore  smootlied  to  tlie  calcu¬ 
lated  spectrum  using  a  moderate  eunount  of  smoothing  in  YOGI.  'I'he  unfolded  and 
calculated  spectra  are  shown  in  Fig.  1.  As  expected  from  the  data,  tlie  SAN4 
matrix  gives  tlie  best  fit.  It  might  bo  argued  that  the  source  encapsulation 
would  broaden  the  spectrum  somewhat,  in  which  case,  the  UTA4  matrix  gives  the 


best  fit.  The  others  show  a  considerable  deficiency  of  high  energy  neutrons. 

The  No.  12  spectrum  is  similarly  compared  with  the  calculated  room- 
return  spectrum  of  Eisenhauer  [12].  A  thermal  component  was  calculated  using 
the  equation  [22] 

nv  =  1.25  Q/S  (5) 

where 

nv  is  the  thermal  neutron  f  luenca  rate  per  cm^ 

Q  is  the  Cf-252  emission  rate  in  neutrons  per  second,  and 
S  is  the  surface  of  the  calibration  room  in  cm^. 

This  component  was  then  combined  with  the  Eisenhauer  calculation  in  order  to 
make  a  comparison  with  the  experimental  data.  Note  in  Table  9  that  the  cal¬ 
culated  counts  from  this  spectrum  also  indicate  that  the  unfolded  spectra  will 
have  lower  energies  when  using  the  M60,  M65,  LOGNM,  and  UTA4  response 
matrices.  Again,  SAN4  and  UTA4  give  good  agreement  with  the  data,  far  better 
than  the  error  of  about  +10%  estimated  for  this  data. 

Using  each  matrix,  spectra  were  unfolded  with  the  calculated 
spectrum  as  a  starting  spectrum,  and  smoothed  to  the  starting  spectrum.  Fits 
were  terminated  when  the  average  error  dropped  below  +1%,  or  after  50 
iterations.  The  spectra  obtained  with  SAN4  and  UTA4  are  shown  in  Fig.  2a. 

The  fit  terminated  after  1  iteration  for  SAN4  and  after  3  iterations  forUTA4. 
As  expected,  UTA4  gives  a  slightly  softer  spectrum.  The  spectrum  obtained 
with  M60,  M65,  and  LOGNM  are  shown  in  Fig.  2b.  Note  the  oscillations  in  the 
M60  and  M65  spectra.  These  are  caused  by  errors  in  the  response  matrices  or 
in  the  detector  data.  LOGNM  gives  a  reasonable  fit  to  the  data;  ±1%  error 
was  reached  after  20  iterations.  However,  the  spectrum  is  considerably  lower 
in  energy  than  the  calculation.  The  data  for  the  BARC  matrix  is  shown  in 
Fig.  2c.  BARC,  M60,  and  M65  did  not  fit  the  data  within  ^1%  after  50  iter¬ 
ations.  The  error  was  approximately  +2%  for  each  matrix.  Using  no  smoottung, 
further  unfolding  produced  oscillations  of  10^  in  the  spectra  without 
significant  improvement  in  the  fit  to  the  data. 
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ig.  1  Cf-252  neutron  spectra  obtained  using  each  of  the 

the  response  matrices  in  the  YOGI  unfolding  code. 

The  spectra  were  unfolded  using  the  NBS  calculation 
as  the  starting  spectrum  for  YOGI. 


LOG  (ENERGY,  MeV) 

Cf-252  room-return  spectrum  obtained  using  the  SAN4 
and  UTA4  response  matrices.  The  calculated 
spectrum  was  used  as  the  starting  spectrum  for  YOGI 
and  the  unfolded  spectra  were  smoothed  to  the 
calculated  spectrum. 
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2b  Cf-252  room-return  spectra  obtained  using  the  M60, 
M65,  and  LOGNM  response  matrices.  The  calculated 
spectrum  was  used  as  the  starting  spectrum  for  YOGI 
and  the  unfolded  spectra  were  smoothed  to  the 
calculated  spectrum. 


LOG  (ENERGY,  MeV) 

Fig.  2c  Cf-252  room-return  spectra  obtained  using  the  BARC 
response  matrix.  The  calculated  spectrum  was  used 
as  the  starting  spectrum  for  YOGI  and  the  unfolded 
spectrum  was  smoothed  to  the  calculated  spectrum. 


Variation  of  Integral  Pareuneters 

Integral  parameters  were  calculated  for  each  spectrum  unfolded  wit^i 
each  matrix  using  Eq.  (3)  and  the  conversion  factors  listed  in  Table  4.  These 
data  are  presented  in  Tcibles  lOa-lOj.  For  each  parameter,  the  average  value 
obtained  with  all  the  matrices  is  given  along  with  the  ratio  of  the  response 
obtained  with  each  matrix  relative  to  the  average  response. 

The  fluence  data  is  shown  in  Table  10a.  This  data  has  been 
normalized  using  spectrum  No.  4,  the  Cf-252  source  with  the  room-return  com¬ 
ponent  removed.  This  simulates  how  a  Bonner  sphere  system  might  typically  be 
calibrated.  Generally,  the  fluence  does  not  vary  more  than  +15%  with  choice 
of  matrix.  The  exception  is  that  M60  tends  to  significantly  overestimate  the 
fluence  for  spectra  having  energies  less  than  the  bare  Cf-252. 

The  data  for  average  energy  is  given  in  Table  10b.  There  is 
reasonable  agreement  for  average  energies  above  approximately  1  MeV.  SAN4  and 
UTA4  give  exceptionally  good  agreement  with  the  average  energies  usually  given 
for  Cf-252  (spectrum  4)  and  AmBe  (spectrum  1)  of  approximately  2.2  and  4.5  MeV. 
Below  this  energy,  there  are  significant  variations  with  SAN4  and  UTA4  giving 
40-70%  over-response  and  LOGNM,  and  especially  BARC,  giving  a  large  under¬ 
response.  The  difference  between  the  average  energy  obtained  with  SAN4  and 
BARC  is  more  than  a  factor  of  10  at  the  lowest  energies.  Spectra  13  through 
20  were  all  generated  with  Cf-252  or  AmBe  in  the  center  of  a  60  cm  dia.  steel 
sphere.  Ing  and  Cross  [23,24]  have  calculated  the  spectrum  for  fission 
neutrons  in  the  center  of  such  a  sphere.  The  calculated  spectrum  gives  the 
same  detector  counts,  within  ±2%,  as  a  Maxwellian  energy  distribution  having 
a  temperature  of  0.30  MeV.  Hence,  we  would  expect  the  Maxwellian  component  of 
spectrum  No.  13  to  have  a  temperature  of  approximately  0.3  MeV.  The  best  fit 
to  the  data  for  SAN4,  M60,  M65,  LOGNM,  BARC,  and  UTA4  for  spectrum  No.  13  are 
0.29,  0.21,  0.27,  0.14,  0.15,  and  0.28  MeV,  respectively.  Adding  Lucite  tends 
to  shift  the  Maxwellian  temperature  upward  [24]  and  t)iis  trend  is  observed. 

We  obtain  a  temperature  of  approximately  0.45  MeV  for  spectra  19  and  20  with 
SAN4  and  UTA4,  with  corresponding  lower  values  for  the  other  matrices. 

Considering  the  variation  in  average  energy,  it  was  expected  that 
there  would  be  considerable  variations  in  dose,  dose  equivalent,  and  quality 
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factor  with  a  strong  correlation  between  the  four  parameters.  This  is  borne 
out  in  Tables  10c,  lOd,  and  lOe  where  the  responses  are  given  respectively  for 
the  dose,  dose  equivalent,  and  quality  factor.  There  are  differences  up  to  a 
factor  of  2  for  dose,  a  factor  of  4  for  dose  equivalent,  and  a  factor  of  2  for 
quality  factor,  with  SAN4  giving  the  highest  values  and  BARC  the  lowest. 

The  calculated  response  per  Rem  of  4  personnel  dosimeters  and  a 
neutron  Remmeter  are  given  in  Tables  lOf-10 j .  For  these  calculations,  the 
response  of  all  devices  was  normalized  to  the  response  obtained  from  the  cal¬ 
culated  Cf-252  spectrum  [21].  The  responses  of  two  albedo  dosimeters  are  given 

in  Tables  lOf  and  lOg.  The  response  of  the  Navy  TLD  badge,  which  consists  of 
7 

a  °LiF-  LiF  thermoluminescence  dosimeter  pair  behind  a  0.38  mm  cadmium  filter, 
is  given  in  Table  lOf.  We  see  significcuit  variations  for  the  lower  energy 
spectra,  with  BARC  giving  significantly  higher  readings  than  the  average,  and 
SAN4  and  UTA4  giving  significantly  lower  readings.  Since  the  conversion 
factors  for  this  table  were  obtained  using  experimental  data  and  the  SAN4 
matrix,  the  actual  response  of  this  badge  to  these  spectra  can  be  obtained  by 
multiplying  the  data  in  column  2  by  the  data  in  column  3.  The  error  on  this 
response  is  approximately  +10%.  Thus,  the  average  calculated  response  for 
these  spectra  varies  by  about  a  factor  of  37,  while  the  measured  variation  is 
about  a  factor  of  20.  This  indicates  that  these  response  matrices  generally 
give  spectra  that  are  too  low  in  energy,  especially  for  spectra  having  low 
average  energies.  The  data  for  a  completely  cadmium  covered  albedo  dosimeter 
is  given  in  Table  lOg.  This  data  is  very  similar  to  that  in  Table  lOF,  however, 
there  is  more  variation  of  the  response  of  this  dosimeter  with  neutron  energy 
than  with  the  Navy  albedo  badge. 

The  response  per  Rem  of  2  personnel  dosimeters  which  have  lower 
energy  cut-offs  is  given  in  Tcibles  lOh  and  lOi.  In  Table  lOh  we  list  the  re¬ 
sponse  of  a  Cr-39  detector  which  has  a  lower  energy  cut-off  of  approximately 
150  keV.  There  are  variations  of  more  than  a  factor  of  6,  depending  on  the 
choice  of  matrix.  The  average  response  of  this  dosimeter  varies  by  about  a 
factor  of  3.3  over  this  energy  range  of  spectra,  with  BARC  showing  the  greatest 
variation,  approximately  a  factor  of  9,  and  M65  the  least,  a  factor  of  approxi¬ 
mately  2.  Our  experimental  data  [25]  gives  a  factor  of  3,  in  good  agreement 
with  SAN4  and  UTA4.  The  response  per  Rem  of  NTA  type  film  is  given  in  Table 


lOi.  Since  NTA  film  has  a  lower  energy  cut-off  of  approximately  0.5  MeV,  we 
might  expect  wide  variations  depending  on  where  the  matrix  tends  to  put  the 
energy  peak  of  the  lower  energy  spectra.  This  is  indeed  the  case;  for  the 
lowest  energy  spectra,  there  is  more  than  3  orders  of  magnitude  difference 
between  BARC  and  the  other  matrices. 

Finally,  in  Table  10 j ,  we  list  the  response  of  the  AN/PDR-70,  an 
Anderson-Braum  type  Remmeter.  The  predicted  response  of  this  instrument  is 
within  ±20%  independent  of  the  choice  of  matrix. 

CONCLUSIONS 

The  choice  of  the  response  matrix  used  to  unfold  neutron  spectra  from 
Bonner  sphere  data  has  a  significant  effect  on  the  spectra  obtained  and  on  the 
integral  parameters  calculated  from  the  spectra.  Our  study  indicates  that  only 
SAN4  and  UTA4  consistently  give  reasonably  shaped  spectra  that  fit  the  sphere 
data  within  experimental  error,  and  agree  well  with  other  calculated  and  ex¬ 
perimental  data.  This  study  indicates  the  need  for  further  calculations  and 
experimental  verification  of  the  Bonner  sphere  response  matrix  in  order  to 
optimize  the  information  derived  from  this  type  of  neutron  spectrometer. 
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TABLE  2 

DESCRIPTIONS  OF  THE  SOURCES  AND  MODERATORS  USED  TO  PRODUCE  THE  NEUTRON  SPECTRA  USED  IN  THIS  STUDY 
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TABLE  3 


THE  DETECTOR  COUNTS,  EXPRESSED  AS  THE  PERCENT  OF  THE  TOTAL  COUNTS, 

FOR  THE  VARIOUS  NEUTRON  SPECTRA 

SPECTRUM 

NUMBER 

DETECTOR  RESPONSE  (%  of  total) 

2  " 

3  " 

5" 

8  " 

10" 

12" 

1 

1 . 04 

5 . 26 

18.70 

28 .01 

25.70 

21 .29 

2 

4 . 25 

9 , 07 

20 . 76 

25.66 

22.26 

18.00 

3 

9 . 10 

13.43 

21 . 64 

22 .78 

18.61 

14 ,44 

4 

1 . 37 

7 . 20 

24.20 

29.13 

22. 52 

15.58 

5 

4 . 89 

12.39 

26.19 

25.15 

18,76 

12.61 

6 

5 .23 

11,59 

25 . 83 

25. 59 

18,93 

12.71 

7 

15.39 

19.34 

23,79 

18.88 

13.55 

9.05 

8 

15.94 

23.00 

26.09 

17.03 

10 . 80 

7,12 

9 

11 .91 

17.65 

25.62 

20,87 

14 .53 

9.42 

10 

16.21 

22.09 

26 . 97 

17.13 

10.87 

6.73 

11 

6,04 

14 . 82 

30,26 

24.45 

15.49 

8.94 

12 

15.91 

24.21 

28 . 80 

16.64 

9 . 30 

5.15 

13 

8 .89 

20.22 

34.29 

20 .73 

10.76 

5.10 

14 

18.66 

27.14 

28 . 87 

14.22 

7 . 54 

3 .57 

15 

11.45 

23.13 

33.61 

18.28 

9.27 

4 . 27 

16 

14. 39 

25 . 38 

32.05 

16.29 

8.27 

3 . 62 

17 

17,53 

26 . 85 

30 . 99 

14.40 

6.92 

3 .31 

18 

17.15 

26 .94 

30 . 64 

14 .80 

7.25 

3 .22 

19 

22,05 

28.93 

28.22 

12.31 

5.85 

2.65 

20 

18.61 

28 . 28 

30.17 

13.47 

6 .49 

2.98 

26 


TABLE  4 

THE  CONVERSION  FACTORS.  RESPONSE  PER  NEUTRON,  FOR  THE 
PARAMETERS  CALCULATED  FROM  THE  UNFOLDED  NEUTRON  SPECTRA 
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MATRIX 


SPECTRUM 


NUMBER 

SAN4 

M60 

M65 

LOGNM 

BARC 

UTA4 

1 

0.60 

2.10 

1.03 

0.68 

7.19 

0.65 

2 

0.63 

5.28 

3.03 

0.50 

6.23 

0.73 

3 

1.35 

1.55 

2.04 

2.55 

5.45 

1.04 

4 

1.97 

1.85 

3.61 

1.58 

11.93 

0.33 

5 

0.57 

1.74 

2.53 

0.78 

7.61 

0.46 

6 

1.28 

2.65 

3.85 

1.01 

7.43 

1.02 

7 

0.52 

2.63 

2.79 

0.70 

5.75 

0.60 

8 

1.31 

2.82 

3.64 

1.47 

3.48 

1.09 

9 

0.96 

3.27 

3.64 

0.49 

6.04 

1.86 

10 

2.13 

3.11 

4.21 

1.35 

4.84 

1.81 

11 

1.30 

4.27 

5.50 

5.44 

8.10 

0.72 

12 

1.23 

3.80 

4.92 

1.39 

4.08 

1.46 

13 

0.90 

3.83 

7.92 

1.82 

8.25 

0.98 

14 

1.18 

3.01 

4.05 

1.57 

4.35 

0.99 

15 

1.03 

3.91 

7.33 

1.56 

6.65 

0.95 

16 

1.47 

3.09 

5.95 

2.29 

5.33 

0.97 

17 

2.24 

5.55 

7.43 

0.89 

3.38 

2.00 

18 

0.71 

3.61 

5.75 

2.01 

2.99 

0.55 

19 

1.18 

4.18 

5.45 

0.93 

1.20 

1.07 

20 

1.09 

4.77 

6.29 

0.57 

2.50 

1.03 

AVERAGE 

1.18 

3.35 

4.55 

1.48 

5.64 

1.02 

S.D. 

0.49 

1.13 

1.88 

1.11 

2.46 

0.46 

Table  6 

The  Average  Error  (%)  on  the  Fit  to  the  Detector  Data 
for  Each  Spectrum  and  Matrix  Using  MAXIET 
and  YOGI  to  Determine  the  Spectrum 


MATRIX 


SPECTRUM 

NUMBER  SANA  M60  M65  LOGNM  BARC  UTAA 


SPECTRUM 

NUMBER 


Table  7a 

The  Percent  Deviation  of  the  Calculated 
Detector  Data  from  the  Experimental 
Detector  Data  for  the  M60  Response  Matrix 


DETECTORS 


0.70 

-1.97 

2.58 

-2.13 

0.65 

0.25 

-0.00 

-0.21 

0.62 

-1.38 

1.34 

-0.35 

-0.19 

-0.17 

0.52 

-1.81 

2,59 

-0.88 

0.65 

-1.26 

1.30 

-2.37 

2.62 

-0.86 

-0.02 

-0.01 

0.27 

-1.65 

2.37 

-0.92 

-0.10 

0.15 

0.17 

-1.83 

2.73 

-1.06 

-0.14 

0.00 

0.53 

-2.07 

2.62 

-0.89 

0.28 

-0.42 

1.28 

-3.56 

3.79 

-1.21 

0.03 

-0.06 

0.66 

-2.50 

2.97 

-1.02 

-0.27 

0.68 

-0.09 

-2.59 

3.84 

-1.45 

-0.03 

0.03 

0.68 

-3.15 

3.97 

-1.37 

0.29 

0,04 

1.37 

-4.56 

4.86 

-1.68 

-0.46 

1.09 

-0.22 

-3.58 

6.60 

-3.11 

-0.02 

0.32 

0.45 

-3.21 

3.72 

-1.13 

-0.08 

0.95 

-0.52 

-3.78 

6.38 

-2.65 

0.05 

0.85 

-0.72 

-2.70 

5.24 

-2.52 

0.37 

2.33 

-2,29 

-3.81 

8.84 

-4.82 

0.14 

0.79 

-0.67 

-3.62 

6.51 

-2.84 

0.45 

1.05 

-1.03 

-3.99 

6.95 

-3.06 

0.35 

1.53 

-1.63 

-3.54 

7.50 

-3.78 

0.10 

0.28 

0.16 

-2.89 

4.30 

-1.77 

0.30 

0.95 

1.12 

0.91 

2.21 

1.27 

The  Percent  Deviation  of  the  Calculated  Detector 
Data  from  the  Experimental  Detector  Data 
for  the  M65  Response  Matrix 


UETIXTORS 


SPECTRUM 

NUMBER  2"  3"  5"  8"  10" 


The  Percent  Deviation  of  the  Calculated 
Detector  Data  from  the  Experimental  Data 
Data  for  the  SAN4  Response  Matrix 


DETECTORS 


SPECTRUM 


NUMBER 

2" 

3" 

5” 

8" 

10" 

12" 

1 

0.22 

-0.63 

0.84 

-0.53 

-0.05 

0.16 

2 

-0.09 

0.37 

-0.63 

0.39 

0.88 

-0.90 

3 

-0.15 

0.28 

-0.32 

-0.05 

0.86 

-0.62 

4 

-0.53 

1.00 

0.22 

-0.95 

-0.25 

0.53 

5 

-0.01 

0.14 

-0.39 

0.73 

-0.59 

0.13 

6 

-0.18 

0.46 

-0.80 

0.80 

-0.03 

-0.24 

7 

-0.16 

0.30 

-0.40 

0.52 

-0.22 

-0.03 

8 

0.22 

-0.49 

0.75 

-1.30 

1.23 

-0.39 

9 

-0.18 

0.39 

-0.61 

0.53 

0.16 

-0.29 

10 

-0.71 

1.53 

-1.86 

1.08 

0.91 

-0.91 

11 

-0.19 

0.48 

-0.53 

0.44 

-0.03 

-0.17 

12 

-0.20 

0.42 

-0.41 

-0.35 

1.22 

-0.67 

13 

-0.15 

0.52 

-0.84 

0.66 

0.09 

-0.27 

14 

-0.04 

0.40 

-0.80 

1.77 

-1.75 

0.46 

15 

-0.17 

0.55 

-0.91 

1.08 

-0.38 

-0.15 

16 

-0.07 

0.20 

-0.25 

1.39 

-2.40 

1.17 

17 

-0.60 

1.67 

-1.78 

0.07 

1.55 

-0.86 

18 

-0.08 

0.24 

-0.37 

0.76 

-0.90 

0.36 

19 

-0.36 

0.99 

-1.25 

0.76 

0.36 

-0.48 

20 

-0.15 

0.60 

-0.85 

0.62 

-0.01 

-0.21 

AVERAGE 

-0.17 

0.47 

-0.56 

0.42 

0.03 

-0.17 

Table  76. 

The  Percent  Deviation  of  the  Calculated  Detector  Data  from 
the  Experimental  Detector  Data  for  the  L06NM  Response  Matrix. 


DETECTORS 


SPECTRUM 


NUMBER 

2" 

3” 

5” 

8” 

10" 

12" 

1 

0.04 

-0.43 

0.92 

-1.03 

-0.05 

0.56 

2 

-0.07 

0.17 

-0.21 

O 

0.75 

-0.49 

3 

0.15 

-0.37 

0.69 

-1.03 

0.49 

0.08 

4 

0.20 

-0.94 

1.97 

-1.44 

-0.68 

0.93 

5 

0.05 

1 

o 

0.43 

0.02 

-0.83 

0.61 

6 

0.09 

-0.27 

0.48 

-0.24 

-0.49 

0.45 

7 

0.10 

-0.22 

0.33 

-0.13 

-0.46 

0.39 

8 

0.34 

-0.63 

1.13 

-1.83 

1.36 

-0.34 

9 

0.14 

-0.25 

0.37 

-0.30 

-0.19 

0.24 

10 

-0.11 

0.27 

-0.33 

0.14 

0.14 

-0.10 

11 

0.31 

-0.79 

1.67 

-1.26 

-0.99 

1.09 

12 

0.17 

-0.43 

0.92 

-1.33 

0.70 

-0.01 

13 

0.15 

-0.58 

1.90 

-2.00 

-1.20 

1.80 

14 

0.36 

-0.86 

1.25 

0.30 

-2.65 

1.66 

15 

0.35 

-1.03 

1.92 

-1.05 

-1.61 

1.47 

16 

0.32 

-0.91 

1.99 

-0.78 

-3.56 

3.07 

17 

-0.11 

0.39 

-0.52 

-0.26 

1.38 

-0.87 

18 

0.06 

1 

o 

1.54 

-1.74 

-2.14 

2.83 

19 

0.10 

-0.50 

1.06 

-1.03 

-0.83 

1.23 

20 

0.09 

-0.35 

0.64 

-0.38 

-0.56 

0.57 

AVERAGE 

0.14 

-0.42 

0.91 

-0.78 

-0.57 

0.76 

S.D. 

0.15 

0.39 

0.78 

0.69 

1.27 

1.03 

Table  7e. 

The  Percent  Deviation  of  the  Calculated 
Detector  Data  from  the  Experimental  Detector  Data 
for  the  BARC  Response  Matrix. 


DETECTORS 


SPECTRUM 


NUMBER 

2" 

3" 

5" 

8" 

10" 

12" 

1 

1.35 

-0.10 

-3.81 

10.45 

-  9.63 

2.87 

2 

-0.10 

2.60 

-4.91 

8.65 

-  8.95 

3.67 

3 

-1.00 

3.55 

-3.50 

5.74 

-  7.74 

3.61 

4 

8.00 

-9.32 

-5.89 

13.43 

-  9.92 

6.19 

5 

0.57 

2.71 

-4.17 

7.79 

-11.24 

5.59 

6 

-0.4  5 

3.90 

-4.88 

8.09 

-10.79 

5.41 

7 

-1.37 

3.10 

-0.62 

3.16 

-  8.28 

4.59 

8 

-0.74 

1.76 

-0.53 

0.56 

-  2.76 

1.79 

9 

-1.31 

3.97 

-1.97 

2.82 

-  8.32 

5.47 

10 

-1.62 

3.79 

-1.57 

1.36 

-  4.98 

3.31 

11 

-0.08 

4.49 

-6.22 

8.02 

-10.40 

5.52 

12 

-1.13 

2.46 

-0.57 

0.15 

-  2.66 

1.83 

13 

-0.53 

7.05 

-4.94 

3.74 

-10.78 

6.74 

14 

-0.93 

1.64 

-0.19 

-0.05 

-  1.75 

1.31 

15 

-1.20 

5.93 

-3.10 

0.62 

-  5.31 

3.49 

16 

-1.68 

4.16 

-0.15 

0.17 

-  4.84 

2.60 

17 

-0.19 

1.13 

-1.75 

0.38 

0.41 

0.04 

18 

-1.04 

2.05 

-0.14 

-0.40 

-  1.98 

1.56 

19 

-0.41 

0.79 

-0.27 

-0.04 

-  0.40 

0.32 

20 

0.31 

-0.11 

-0.99 

0.22 

0.85 

-0.26 

AVERAGE 

-0.18 

2.28 

-2.51 

3.74 

-  5.97 

3.28 

S.D. 

2.07 

3.29 

2.12 

4.25 

4.11 

2.12 

34 


Table  7f 


» 

The  Percent  Deviation  of  the  Calculated  Detector  Data 
from  the  Experimental  Detector  Data  for  the  UTA4  Response  Matrix. 


DETECTORS 
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Table  9. 

The  Measured  Counts,  Calculated  Counts,  and  Percent  Difference  Between 
Calculated  Counts  and  Measured  Counts  for  the  Room-Return  Spectrum  of  Cf-252. 


Table  10a. 

The  Average  Neutron  Fluence  Obtained  Using  All  the  Matrices 
and  the  Ratio  of  the  Fluence  Obtained  with  each 
Matrix  to  the  Average  Fluence 
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